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We report the result of our ab initio calculation of the &s 2 S\/2 — * 5d 2 Ds/2 parity nonconserving 
electric dipole transition amplitude in 137 Ba + based on relativistic coupled-cluster theory. Consid- 
ering single, double and partial triple excitations, we have achieved an accuracy of less than one 
percent. If the accuracy of our calculation can be matched by the proposed parity nonconservation 
experiment in Ba + for the above transition, then the combination of the two results would provide 
an independent non accelerator test of the Standard Model of particle physics. 



Parity nonconservation (PNC) in atoms arising from 
neutral weak currents has the potential to test the Stan- 
dard Model (SM) of particle physics 0,0 By combining 
the results of high precision measurements and calcula- 
tions of atomic PNC observables, it is possible to ex- 
tract the nuclear weak charge Q and compare with its 
corresponding value in the SM. A discrepancy between 
these two values could reveal the possible existence of 
new physics beyond the SM. The most accurate data on 
atomic PNC currently comes from the Qs 2 S\/2 — > 7s 2 Si/2 
transition in cesium (Cs), where the claimed experimen- 
tal Q and theoretical Q accuracies are 0.35% and 0.5% 
respectively and the deviation from the SM is about la 
4]. It would indeed be desirable to consider other can- 
didates which could yield accurate values of the nuclear 
weak charge. In this context an experiment to observe 
PNC in the Qs 2 S\/2 — > 5d 2 D 3 / 2 transition in Ba + using 
the techniques of ion trapping and laser cooling proposed 
by Fortson is of special importance pa, Ig • 

This Letter is concerned with a high precision calcu- 
lation of the amplitude of the above mentioned parity 
nonconserving electric dipole (Elpjvc) transition in Ba + 
using relativistic coupled-cluster (RCC) theory, which is 
equivalent to all order relativistic many-body perturba- 
tion theory 0- It is the first application of this the- 
ory to atomic PNC. Blundell et al had used this theory 
in the linear approximation to calculate Elpjvp for the 
6s 2 S]/ 2 — > 7s 2 Si/ 2 transition in Cs 8J. Dzuba et al 
and Geetha have calculated this PNC amplitude for 
the 6s 2 5i/2 — * 5d 2 D 3 /2 transition in Ba + as discussed 
later. 

The parity nonconserving nuclear spin independent 
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(NSI) interaction arises from the nucleon-electron neu- 
tral weak interaction and its Hamiltonian is given by 

H ™ c = £k Qw 7sMr) ((U) 

where Gf is the Fermi constant, Q w is the nuclear weak 
charge which is equal to [(2Z + N)ci u + (27V + Z)c\J[ with 
c\ u and cid representing electron-up-quark and electron- 
down-quark coupling constants respectively, pn{t) is the 
nuclear density function and Js(= 270717273), which is a 
pseudo-scalar, is the product of the four Dirac matrices. 

Hpnc is responsible for mixing atomic states of oppo- 
site parities but with the same angular momentum. Its 
strength is sufficiently weak for it to be considered as a 
first-order perturbation. It is therefore possible to write 
the n'th state atomic wavefunction as 

\9 n ) = \9W)+G F \*W). (0.2) 

In RCC, the atomic wavefunction |\E , !°'') for a single va- 
lence (v) open-shell system is given by fRUlf 

|^)) = e T<0, {l + ^}|$„) (0.3) 

where we define |$„) = aJ|$o)i with |$q) as the Dirac- 
Fock (DF) state for closed-shell system. 

In the singles and doubles approximation we have 

T (o) _ T (o) + T (o) 

= s[J + S 2 J (0-4) 

where t[ 0) and T 2 (0) are the single and double particle- 
hole excitation operators for core electrons and and 
S$ are the single and double excitation operators for 
the valence electron respectively. The amplitudes corre- 
sponding to these operators can be determined by solving 
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the relativistic coupled-cluster singles and doubles equa- 
tions. A subset of important triple excitations have been 
considered in the determination of the open shell ampli- 
tudes S$ and which is described in 0,0. 

Using eqn. (0.2), the explicit form of Filp^c, is given 

by 



El 



PNC 
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(0.5) 



where D is the electric dipole (El) operator, i and / 
subscripts are used for initial and final valence electrons 
respectively. Using the explicit expression for the first 
order perturbed wavefunction, we get 



(^°)| D |^)(*ri H PNci* 

Ei — Ej 

(^ 0) |HNSi c |^) )( ^)|p|^ 
Ef-Ei 



(0)\/, t ,(0)|ttNSI h T /(°)\ 



■E 



(0.6) 



where / represent intermediate states. 

It is obvious from the above equation that, the accu- 
racy of the calculation of Elpjvc depends on the exci- 
tation energies of the different intermediate states, the 
matrix elements of Hpj^ c and D. Blundell et al have 
used the above equation to determine ElpNc for the 
6s 2 Si/2 — * 7s 2 Si/2 transition in Cs by considering the 
most important intermediate states @ . The drawback of 
this approach is that the summation can be performed 
only over a finite set of intermediate states which limits 
the accuracy of the calculation. The method we have 
used in the present work circumvents this problem by 



solving the first order perturbed equation 

( H (o) _ £(0)^(1)) = (£(D „ H^ c )|¥<°>>. (0.7) 

where E^ vanishes for first order correction. 

The perturbed cluster operators can be written as 



T = + G f TW, 
S v = S { v 0) + GpSP 



(0.8) 



where and Sv are the first order Gp corrections to 
the cluster operators and S^ ' respectively. The am- 
plitudes of these operators are solved, keeping up to liner 
in PNC perturbed amplitudes, by the following equations 



(^IH^TW+HNSI |$ ) = 0, 



N 



H»$o> 



0, 



(0.9) 



and 



(*5|H$S« + (i4 0) T« + H nsi c){ i + £C°)}|$ fl ) 

= -<*g|S-W|$„)IP, 



(Kll^s^ + (h^tw + hnsi c){ i + S (0) } ^ v) 

= -«I|S«|^)IP(0.10) 

where IT- ' is the Dirac-Coulomb (DC) Hamiltonian and 
H is defined as e _T< ° He T<0> which is computed after 
determining T^°\ IP is the ionization potential energy 
corresponding to the valence electron V and the sub- 
script N denotes normal form of an operator. We have 
used a, b.. and p, q.. etc. to represent holes and particles 
respectively. and |$^) are the single and double 

excited states respectively with respect to |3>«). Using 
Eqns. (0.3), (0.5), (0.8) and only keeping terms linear in 
Gf, the expression for Elpjvc can be written as 



< * f \{l + sfV +T( 1 )V (0)t De T<0) {l + T( 1 ) + T^S^ +5f ) }|* i > 

ElpTVC = , = 

^/(1 + Nf)(l + N^) 

< $ f \S?>' B®(1 + S| 0) ) + (1 + 5|. 0)t )DW5f 5 + sf'^D'DM+pTWjsf 1 + (Tf^'DW + PrWjsf^, > 
- = 

+ )(l + ivf 0) ) 



In the above expression we define D(°) = e T De T and N{, = Sv e T e T Si for the valence electron V and 
each term is connected. The above matrix element is evaluated by a method similar to that used in our earlier works 
of Ba+ mEt. 



The orbitals are constructed as linear combinations of 
Gaussian type orbitals (GTOs) of the form 0] 

Fi, k (r) =r k e- a ^. (0.12) 
where k — 0,1,.. for s,p,.. type orbital symmetries re- 



spectively. For the exponents, we have used 

on = (0.13) 

We have considered 30si/2, 25pi/2, 25p 3 / 2 , 25d 3 / 2 , 
25d 5/2 , 20/ 5/2 , 20/ 7/2 , 20 ffr/2 and 20 59/2 GTOs for the 
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TABLE I: Excitation energy (cm -1 ), El transition ampli- 
tudes (a.u.) and magnetic dipole hyperfine structure constant 
(MHz) for different low-lying states of Ba + . 



TABLE II: Square root of the magnetic dipole hyperfine con- 
stants (MHz) and their deviations from experimental results. 

Experiment This work Deviation (%) 



Initial state 
-^Final state 



Excitation 
energy 
Expt. [17] 



El transition 
amplitude 
Expt. [18] 



Atomic state 



Hyperfine 
constant (A) 
Expt. [19, 20, £1] 



6s 2 5* 1/2 
6p 2 Pi/ 2 



6s 2 S 1/2 
6/P3/2 



5d 2 D 



3/2 



5d 2 D 



3/2 



QP*PU2 6p"P 3 /2 



20410 22104 15097 16795 
20262 21952 15388 17079 



\ J A 6s 2 S 1/2 A 6y 2 P 1/2 
\/ A 6p 2 P 3/2 A Sd 2 D 3/2 



1728.83 
155.35 



1738.1 
156.08 



0.5 
0.5 



3.37 4.72 3.08 1.36 

3.36(0.16) 4.67(0.08) 3.03(0.08) 1.36(0.04) 



TABLE III: Contributions 
xlO _11 jeao(— Qw/N) using 



to the El pnc calculation in 
RCC calculation. 



6s 2 S 1/2 6p 2 P 1/2 6p 2 P 3/2 5d 2 D 3 /2 



Initial pert, 
terms 



4078.18 740.77 128.27 189.92 
4018.871(2) 743.7(3) 127.2(2) 189.7288(6) 



Qs 2 S[ L ; 2 - 
5d 2 D^ 



Final pert, 
terms 



Dirac-Fock contribution 



6s 2 S[ u > 2 - 
5rf 2 ^/2 



DH 



NSI 
PNC 



2.018 



-0.3 xlO" 



DF calculation and all occupied (active holes) orbitals 
in the RCC calculations. We have chosen ao as 0.00525 
and (3 as 2.73 for all the symmetries. All orbitals are 
generated on a grid using a two-parameter Fermi nuclear 
distribution approximation given by 



P = 



Pa 



1 -f e {r- c )/ a 



(0.14) 



where po is the average nuclear density, 'c' is the half- 
charge radius, and 'a' is related to the skin thickness. 

Our earlier calculations of excitation energies 0], El 
transition amplitudes and magnetic dipole hyper- 
fine constants 0] for some of the low-lying states in 
Ba + based on RCC theory suggest that it is in princi- 
ple possible to perform a calculation of Elpjvc for the 



6s 2 S 1 



/2 



5d -D3/2 transition in that ion to an accu- 



racy of better than one percent. We have recalculated 
these quantities using the same method but with a larger 
basis and the results are given in table I. The agree- 
ment with experiment of the most important excitation 
energy (6p 2 Pi/ 2 ) for the calculation of EIpnc is less 
than one percent. This is also the case for the hyper- 
fine constants of three of the states - Qp 2 Pi/ 2 , 6p 2 P 3 / 2 
and 5d 2 D 3 / 2 , while for the 6s 2 Si/ 2 state, the agreement 
is a little over one percent. All the transition ampli- 
tudes are within the experimental error bars. The result 
of our calculation of the electric quadrupole (E2) ampli- 
tude for the 6s 2 Si/2 — > 5d 2 D 3 / 2 transition is 12.61 in 
a.u. It is in agreement with our earlier calculation [T^ | 
and well within the experimental bounds p^ . In ta- 
ble II, we present the values of the square root of the 
product of the hyperfine constants. The accuracies of 
these two quantities give an indication of the accuracies 
of the PNC matrix elements between 6s 2 S' 1 / 2 and 6p 2 P\/ 2 
states as well as 6p 2 P 3 / 2 and 5d 2 D 3 / 2 states. Both of 
them are in excellent agreement with experiment, sug- 
gesting that the two leading PNC matrix elements used 
in the Elp^vc calculation are very accurate. The contri- 
butions from the different terms in EIpatc" are presented 
in table III. It is clear that the largest contribution comes 



DOTS™ 



(1) 



Norm. 



Total 



0.0003 
2.634 
-0.242 
0.149 
0.007 
-0.116 
-0.001 
-0.046 

2.375 



rW T D 



(0) 



s^ 3t b(vs\ 



(0) 
li 
(0) 



0.418 
-0.179 
-0.166 
0.003 
0.008 
-0.009 
0.001 
-0.001 

0.087 



from DS' 1 1 ' > (diagram l(m)) which represents the DF term 
and a certain sub class of core polarization as well as 
pair correlation effects This is due to the relatively 



large (6si/ 2 — 6pi/ 2 ) cluster amplitude. Two dif- 
ferent types of core polarization effects; DT^ (diagram 
l(i)) and DS^ 1 ' as well as its conjugate (diagrams l(v) 
and l(vi)) also make significant contributions. The former 
is mediated by the neutral weak interaction and involves 
the 6s valence and core electrons. Correlation effects cor- 

and S^^DS^ are non negli- 



responding to S^^DS^ 1 ' 
gible, but their signs are opposite. Contributions from 
other terms are comparatively small. 

The result of F,1 PNC for the 6s 2 5i/ 2 — > hd 2 D 3 / 2 tran- 
sition in our calculation is 2.46 x 10~ n ieao(— Qw/N). 
It is larger in magnitude than those obtained by Dzuba 
et al and Geetha 01 as shown by table IV. The for- 
mer work is based on a variant of all order many-body 
perturbation theory, but it has some semi-empirical fea- 
tures. It is carried out by using two different approaches. 
One of them is similar to the sum-over-states approach 
by Blundell et al Q and the other is known as the mixed 
approach where the PNC interaction explicitly mixes 
states of opposite parities. However, both calculations 
do not include contributions from certain correlation ef- 
fects; i.e. structural radiation, weak correlation potential 
and normalization of states lij that are included in our 



calculation. Their 6p P 1 



5d D 3 / 2 El matrix ele- 
ment which is important for the above mentioned PNC 



/2 
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FIG. 1: Important Goldstone diagrams corresponding PNC 
amplitudes. 




(v) <vi) (vii) (viii) 



TABLE IV: Comparison of EIpmc results from different cal- 
culations in xlO~ n ieao(— Qw/N). 



Dzuba et al [9] 
(mixed parity) (sum-over-states) 


Geetha [10] 


Present work 


2.17 2.34 


2.35 


2.46 ± 0.02 



transition amplitude is not as accurate as ours. Fur- 
thermore, the accuracies of their PNC matrix elements 
are not known as they have not performed calculations 
of the hyperfine constants of the relevant states. The 
reason for the discrepancy between our calculation and 
Geetha's is that our approach implicitly includes several 
intermediate states; particularly doubly excited opposite 
parity states which her sum-over-states approach omits. 



The error accrued in our calculation of Elp^c can be 
determined from the errors in the excitation energies, El 
transition amplitudes and hyperfine constants (see table 
I). We have not estimated the errors in the calculated val- 
ues of these quantities by comparing with measurements, 
since the error bars in the El transition amplitudes are 
rather large. Instead, we have taken the differences of our 
RCC calculations with single, double and leading triple 
excitations and just single and double excitations as the 
errors. The error in Elpjvc (0.02) has been obtained by 
adding the errors for the different quantities it depends on 
in quadrature for the leading intermediate states &p 2 P\/2 
and 6p 2 P 3 / 2 and using a scale factor to estimate the er- 
rors from other intermediate states that together make a 
small contribution. 

The contribution of the Breit interaction to El pnc a t 
the DF level is 0.1% and the nuclear structure contri- 
bution is 0.3%. The latter has been determined more ac- 
curately than Blundell et al. Q using relativistic mean 
field theory. 

In conclusion, we have performed a sub one percent cal- 
culation of Elpjvc for the 6s 2 Si/ 2 — > 5d 2 D 3 / 2 transition 
in Ba + using RCC. We have included single, double as 
well as a leading class of triple excitations and highlighted 
the importance of various many-body effects. Given the 
promise that the Ba + PNC experiment holds out, it does 
indeed appear that in the future the result of that exper- 
iment combined with our calculation would constitute a 
new and an important probe of physics beyond the SM. 
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